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. Abstract 

, Integrable models are often constructed with real systems in mind. The exact solvability of the 
models leads to results which are unambiguous and provide the correct physical picture. In this review, 

^-H i we discuss the physical basis of some integrable spin models and their relevance in the study of real 

^ . systems. The emphasis in the review is on physical understanding rather than on the mathematical 

0^ ' aspects of integrability. 

m : 
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22 ' 1 Introduction 
O ; 

' The study of integrable models constitutes an important area of theoretical physics. Integrable models in 

I condensed matter physics describe interacting many particle systems. The most prominent examples are 

' interacting spin and electron systems which include several real materials of interest. Integrable models, 

^ I because of their exact solvability, provide a complete and unambiguous understanding of the variety of 

I ' phenomena exhibited by real systems. Integrability in the quantum case implies the existence of N conserved 

' ^ I quantities where N is the number of degrees of freedom of the system. The corresponding operators including 
the Hamiltonian commute with each other. More specifically, integrable models are also described as exactly- 



Q [ solvable since the ground state energy and the excitation spectrum of the models can be determined exactly. 

Historically, the first example of the exact solvability of a many body problem was that of a spin— ^ quantum 
spin chain ||l||. The technique used to solve the eigenvalue problem is now known as the Bethe Ansatz (BA) 
' named after Hans Bethe who formulated it. The demonstration of integrability, namely, the existence of N 
H I commuting operators can be made in the more general mathematical framework of the Quantum Inverse 
Scattering Method (QISM) The BA has been used extensively to obtain exact results for several quantum 
models in one dimension (Id). Examples include the Fermi and Bose gas models in which particles on a 
line interact through delta function potentials |Q, the Hubbard model |Q, Id plasma which crystallises as a 
Wigner solid |^, the Lai-Sutherland model which includes the Hubbard model and a dilute magnetic model 
as special cases |^, the Kondo model the single impurity Anderson model |^, the supersymmetric t — J 
model (J — 2t ) etc The BA method has further been applied to derive exact results for classical lattice 
statistical models in 2d. 

The BA denotes a particular form for the many-particle wave function. In a Id system with pairwise 
interactions, a two particle scattering conserves the momenta individually due to the energy and momen- 
tum conservation constraints peculiar to Id. Hence the scattering particles can either retain their original 
momenta or exchange them. In the case of two particles {N = 2), the wave function has the form 

where xi , X2 denote the locations of the two particles and fci , are the momentum variables. The wave 
function can alternatively be written as 

i!{xi,X2) = e^('=i^i+''^^^) + e*»i2e*('=2^i+'=i^2) (-2) 
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where 6*12 is the scattering phase shift. The BA generahses the wave function (Eq.(l)) to the general case of 
N particles and is given by 

EN 
^=1 <2;2 < xn (3) 

p 

The sum over P is a sum over all permutations of 1, TV. The amplitude A{P) is factorisable. Each A{P) is 
a product of factors e^'J 's corresponding to each exchange of ki 's required to go from the ordering 1, iV to 
the ordering P. An overall sign factor may arise depending on the parity of the permutation. The unknown 
variables 0ij 's and ki 's are obtained as solutions of coupled nonlinear equations. The factorisability condition 
is at the heart of the exact solvability of the eigenvalue problem. In the more general QISM approach, the 
so-called Yang-Baxter equation provides the condition for factorization of a multi-particle scattering matrix 
in terms of two-particle scattering matrices. 

The traditional BA (Eq.(3)) is known as the Coordinate Bethe Ansatz (CBA). Over the years, the BA 
method has been generahsed in different ways. The nested BA technique Q has been applied to study a 
system of particles with internal degrees of freedom. The state of a system of electrons is specified in terms 
of both the spatial positions as well as the spin indices of the electrons. The Asymptotic Bethe Ansatz 
deals with a class of models in which the interaction between a pair of particles falls off as the inverse square 
of the distance between the particles. The Thermodynamic Bethe Ansatz method (l^ is used to calculate 
thermodynamic quantities and is a finite temperature extension of the BA method. The Algebraic Bethe 
Ansatz (ABA) has been developed in the powerful mathematical framework of the QISM. The ABA and 
CBA are equivalent in the sense that both lead to the same results for the energy eigenvalues. The CBA, 
however, does not provide knowledge of the correlation functions as the structure of the wave function is 
not sufficiently expHcitly known. The QISM allows the calculation of the correlation functions in some cases 
. The mathematical formalism is also much more systematic and general. One can further estabHsh the 
existence of an infinite number (TV 00) of mutually commuting operators. The QISM moreover provides a 
prescription for the construction of integrable models. In this review, we will not discuss the mathematical 
aspects of integrable models for which a good number of reviews already exist |§ ||, 0- We focus 
on the physical basis of some integrable spin models in condensed matter physics and the useful physical 
insights derived from the solution of these models. The review is not meant to be exhaustive and should be 
supplemented by the references quoted at the end. 

2 Spin models in Id 

The interest in Id spin models arises from the fact that there are several real magnetic materials which 
can be described by such models. The spins interact via the Heisenberg exchange interaction and in many 
compounds the exchange interaction within a chain of spins is much stronger than that between chains. 
Thus the compounds effectively behave as linear chain systems. The most general exchange interaction 
Hamiltonian describing a chain of spins in which only nearest-neighbour (n.n.) spins interact is given by 

N 

HxYZ = [-^-^rsr+i + JySlSl^, + J.StSt^^ (4) 

where Sf{a — x,y,z) is the spin operator at the lattice site i, N is the total number of sites and Ja 
denotes the strength of the exchange interaction. Consider the spins to be of magnitude ^. The eigenvalue 
problems corresponding to the isotropic chain {J^ — Jy = Jz = J) and the longitudinally anisotropic chain 
{Jx = Jy 7^ Jz) were originally solved using the CBA. Later, the same solutions were obtained using the 
formahsm of QISM [|l3[ Baxter [^calculated the ground state energy of the fully anisotropic model 
(Eq.(4)) and Johnson, Krinsky and McCoy |[l9j found the excitation spectrum. The results were derived 
on the basis of a special relationship between the transfer matrix of the exactly-solved 2d classical lattice 
statistical eight vertex model and the fully anisotropic quantum spin Hamiltonian Hxyz ■ Later, the same 
results were obtained by the ABA approach of the QISM. The Ising {Jx — Jy = 0) and the XY {Jz = 0) 
Hamiltonians are special cases of Hxyz ■ 
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Consider the isotropic Heisenberg exchange interaction Hamiltonian in Id 



N 

H = JY,^^■'^^+l (5) 

i=l 

with periodic boundary conditions. The sign of the exchange interaction determines the favourable ahgnment 
of the n.n. spins. J > corresponds to antiferromagnetic (AFM) exchange interaction due to which n.n. 
spins tend to be antiparallel. If J < (equivalently, replace J by — J in Eq.(5) with J > ), the exchange 
interaction is ferromagnetic (FM) favouring a parallel ahgnment of n.n. spins. One can include a magnetic 
field term —h X^^Li in the Hamiltonian (Eq.(5)), where h is the strength of the field. Given a Hamiltonian , 
the quantities of interest are the ground state energy and the low- lying excitation spectrum. Knowledge of the 
latter enables one to calculate thermodynamic quantities like magnetization, specific heat and susceptibility 
at low temperatures. In the case of the FM Heisenberg Hamiltonian, the exact ground state has a simple 
structure. All the spins are parallel, i.e., they align in the same direction. The lowest excitation is a spin 
wave or magnon. The excitation is created by deviating a spin from its ground state arrangement and 
letting it propagate. For more than one spin deviation, one has continua of scattering states as well as 
bound complexes of magnons. In a bound complex, the spin deviations preferentially occupy n.n. lattice 
positions. The r-magnon bound state energy can be calculated using the BA and the energy (in units of 
J) measured w.r.t. the ground state energy is 

e=-{l-cosK) (6) 
r 

where K is the centre of mass momentum of the r magnons. The spin wave excitation energy is obtained for 
r = 1. The results can be generahsed to the longitudinally anisotropic XX Z Hamiltonian. The multimagnon 
bound states were first detected in the quasi-ld magnetic system CoCl2-'2.H20 [|o|. Later improvements 
made it possible to observe even 14 magnon bound states pi]| . 

In the case of the AFM isotropic Heisenberg Hamiltonian, the ground state is a singlet and the ground 
state wave function is a linear combination of all possible states in which half the spins are up and the 
other half down. The AFM ground state can be obtained from the FM ground state by creating r = ^ 
magnons with momenta ki and negative energies — J(l — cosh). Remember that the sign of the exchange 
integral is changed in going from ferromagnetism to antiferromagnetism. The highest energy state in the 
FM case {r = y) becomes the ground state in the AFM case. The BA equations can be recast in terms of 
the variables Zi = cot{^-) ||2^ : 

Narctanzi = nli + arctan ^ ^ 2 ) ' * ^ ^' ^ ('^) 

The Bethe quantum numbers /i's are integers (half integers) for odd (even) r. For a state specified by 
{/i, /r}, the solution {zi,...,Zr) can be obtained from Eq. (7). The energy and the momentum wave 
number of the state are given by 



J 



i=l 



i=l 

with Ep = For the AFM ground state, the Bethe quantum numbers are given by 

In the thermodynamic limit N ^ 00, the exact ground state energy has been computed as 

1 



Eg = NJi-ln2+-) (11) 
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The AFM ground state serves as the physical vacuum for the creation of elementary excitations. These 
excitations are not the spin-1 magnons but spin— ^ spinous ||2^. The spinons can be generated systematically 
by suitable modifications of the vacuum array of the BA quantum numbers (Eq.(2)) (for details see ||2^, ^ ). 
For even N, spinons are always created in pairs, each such pair originating from the removal of one magnon 
from the ground state. Since the spinons are spin— ^ objects, the lowest excitations consisting of a pair of 
spinons are four-fold degenerate , three triplet {S = 1) and one singlet {S = 0) excitations. The energy can 
be written as E{ki, = e(fci) + £(^2) where the spinon spectrum e(fci) = § sinfe^ and the total momentum 
k = ki + k2. At a fixed total momentum k, one gets a continuum of scattering states. The lower boundary of 
the continuum is given by ^ \sink\ with one of the fc-s = 0. The upper boundary is obtained for fci = fc2 = | 
and is given by tt |sm||. Figure 1 gives an example of a two-spinon configuration. 



Figure 1. A two-spinon configuration in an AFM chain. 

The BA results are obtained in the thermodynamic limit. In this limit, the energies and the momenta 
of the spinons just add up, showing that they do not interact. Since the spinons are excited in pairs, the 
total spin of the excited state is an integer. Inelastic neutron scattering study of the linear chain S = ^ 
HAFM compound KCuF^ has confirmed the existence of unbound spinon pair excitations ||2^. It is to be 
noted that in the case of a ferromagnet, the low-lying excitation spectrum consists of a single magnon branch 
whereas the AFM spectrum is a two-spinon continuum with well-defined lower and upper boundaries. 

The dynamical properties of a magnetic system are governed by the time-dependent pair correlation 
functions or their space-time double Fourier transforms known as dynamical correlation functions. An 
important time-dependent correlation function is 

G(i?,t) = (^fl,(0- 5^0(0)) (12) 

The corresponding dynamical correlation function is the quantity measured in inelastic neutron scattering 
experiments. The differential scattering cross-section in such an experiment is given by 

^2 1 _y r+oc 

^ (X S'^^it^u;) ^-Y^e^qR j ^ dte^-' (^^^W^o^W) (13) 

where if and uj are the momentum wave vector and energy of the spin excitation and /i — x,y,z. For a 
particular 'q' , the peak in S^^^{lf,uj) occurs at a value of uj which gives the excitation energy. At T = 0, 

S'^'^ir, ^)^Y1 '^(^ + Ex) (14) 

A 

Eg{E\) is the energy of the ground (excited) state and 

Mj: = 27r|(G|5^(Y)|A)|^ (15) 

is the transition rate between the singlet (5404 — 0) ground state \G) and the triplet {Stot — 1) states |A) | [2^ . 
The exact calculation of the dynamical correlation functions in the BA formaHsm is not possible. Bougourzi 
et al have used an alternative approach, based on the algebraic analysis of the completely integrable 
spin chain, and have calculated the exact 2-spinon part of the dynamical correlation function S^^((i,ui) for 
the Id 5 = 5 AFM XXZ model. In this model, the Ising part of the XXZ Hamiltonian provides the 
dominant interaction. Karbach et al have ||2^calculated the exact 2-spinon part of S^^{q,iL!) for the isotropic 
Heisenberg Hamiltonian. In both the cases, the size of the chain is infinite. The exact form of the 2-spinon 
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contribution to the dynamical correlation function S^^{q,uj) of the S 
and is given by 



XX Z HAFM chain is complicated 










-^0 



(16) 



where / = ^sinh^^' - - i-*^ 
K(k'). The anisotropy parameter g 



X ' X = f+p"cind fc, /c' = Vl— A? are the modulii of the elliptic integrals K = K{k), K' = 



-exp{- 



-■kK' 
K 



) with -i = A = (2±f— ). Also, 



2/sm((7) 



/31('Z,'^) 

(F is the incomplete elliptic integral) 



Wo 

1 + X 



1 



«^cs«n — — r-^ ,X 



2 , e'^'- cosh(2jl)cos(t^l) - 1 



(=1 



I sinh{2^l)cosh(-^l) 



(17) 
(18) 
(19) 

(20) 
(21) 



7 = ^^,t = and z9ti(a;) is a Neville theta function. The derivation of S^2){1t^) involves generat- 
ing the 2-spinon states from the spinon vacuum, namely, the AFM ground state, with the help of spinon 
creation operators and expressing the spin fluctuation operator S^{q) in terms of the spinon creation op- 
erators. The 2-spinon part is expected to provide the dominant contribution to the dynamical correla- 
tion function (Eq. (14)). For example, in the case of the isotropic Heisenberg Hamiltonian, the 2-spinon 
excitations account for approximately 73% of the total intensity in S''^(q,uj). The 2-spinon triplet ex- 
citations play a significant role in the low-temperature spin dynamics of quasi-ld AFM compounds like 
KCuF3,Cu{C6D5COO)2.SD20,Cs2CuCh and Cu(C4-ff4A^2(A^03)2) H |2|. These excitations can be 
probed via inelastic neutron scattering and hence a knowledge of the exact dynamical correlation function 
is useful. The 2-spinon singlet excitations cannot be excited in neutron scattering because of selection rules 
(the spinon vacuum \G) is a singlet and the excited state |A) in Eq.(15) is a triplet). Linear chain compounds 
like CuGeOs exhibit the spin-Peierls transition ||2^. The transition gives rise to lattice distortion and con- 
sequently to a dimerization of the exchange interaction. Exchange interactions between successive pairs of 
spins alternate in strength. There is a tendency for the formation of dimers (singlets) across the strong 
bonds. One can construct an appropriate dynamical correlation function in which the dimer fluctuation 
operator (DFO) replaces the spin fluctuation operator S^{q). The DFO connects the AFM ground state to 
the 2-spinon singlet and not to the 2-spinon triplet. 

Two well-known physical realizations of the Id = ^ Ising-Heisenberg compounds are CsCoCl^ and 
CsCoBr^ . Several inelastic neutron scattering measurements have been carried out on these compounds to 
probe the low-temperature spin dynamics ||3^. In these compounds, the Ising part of the XX Z Hamiltonian 
is signiflcantly dominant so that perturbation calculations around the Ising limit are feasible. Near the Ising 
limit, the exact 2-spinon dynamical correlation function S^'^{q,uj) is identical in the lowest order to the 
first-order perturbation result of Ishimura and Shiba (IS) . The IS calculation provides physical insight 
on the nature of spinous. The Ising part of the XXZ Hamiltonian is the unperturbed Hamiltonian and 
the XY part constitutes the perturbation. The two-fold degenerate Neel states are the ground states of the 
Ising Hamiltonian. These two states serve as the "spinon vacuua". An excitation is created by flipping a 
block of adjacent spins from the spin arrangement in the Neel state. For example, in Figure 1, a block of 
seven spins is flipped in the Neel state. The block of overturned spins gives rise to two parallel spin pairs 
at its boundary with the unperturbed Neel configuration. It is these domain walls or kink solitons which 
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are the equivalents of spinons. A 2-spinon excited state {S^^^ = 1) is obtained as a linear superposition of 
states in which an odd number v {v = 1,3,5, ....) of spins is overturned in the Neel configuration. In each 
such state, both the domain walls have equal spin orientations with the spins pointing up. The excitation 
continuum of two spinons is obtained in first order perturbation theory. The Hneshapes of S^^{q, uj) observed 
in experiments are highly asymmetric with a greater concentration of intensity near the spectral threshold 
and a tail extending to the upper boundary of the continuum. The exact 2-spinon part of S^^{q, w) has also 
an asymmetric shape in agreement with experimental data. The first order perturbation-theoretic result of 
IS for S^^{q,uj) fails to reproduce the asymmetry. A second-order perturbation calculation leads to greater 
asymmetry in the lineshapes Furthermore, in the framework of a first order perturbation theory, the 
effects of full anisotropy {Jx ^ Jy ^ Jz), next-nearest-neighbour coupHng, interchain coupling and exchange 
mixing have been shown to give rise to asymmetry in Hneshapes fs^ . 

Recently, a large number of studies have been carried out on a class of models in which the interaction 
between spins falls off as the inverse-square of the distance between them. A lattice model which belongs to 
this class is known as the Haldane-Shastry model [|4j the Hamiltonian of which is given by 

where d{l) = (^) |sm^|is the chord distance between the pair of spins separated by I sites on a ring with N 
equally spaced spins. Pij is the spin exchange operator, Pij — {2~^i.~^ j + ^). The model is exactly solvable 
and the key results are: the ground state has a form similar to the fractional quantum Hall ground state, the 
ground state is a QSL and the elementary excitations are the spin— ^ spinons obeying fractional statistics, 
the thermodynamics as well as the various dynamical correlation functions can be calculated exactly. The 
latter calculations are possible because of the simple structure of the eigenspectrum. 

A correct analysis of the BA equations for the S = ^ HAFM in Id gave rise to the concept of spinons which 
has subsequently been verified in experiments. Approximate methods like spin wave theory fail to predict 
the spinon continuum thus pointing to the importance of integrable models in providing the correct physical 
picture. The existence of spinons in dimension greater than one is a highly debatable issue. No precise 
statement can be made due to the lack of exact results in c? > 1. The issue is of considerable significance 
in connection with the resonating-valence-bond (RVB) theory of high temperature superconductivity. In a 
valence bond (VB) state, pairs of spins are in singlet spin configurations (a singlet is often termed as a VB). 
The RVB state is a coherent Hnear superposition of VB states. In 1973, Anderson |Q in a classic paper 
suggested that the ground state of the S — ^ HAFM on the frustrated triangular lattice is a RVB state. 
The RVB state is a singlet (total spin is zero) and is often described as a quantum spin liquid (QSL) since 
translational as well as rotational symmetries are preserved in the state. The RVB state is spin disordered and 
the two-spin correlation function has an exponential decay as a function of the distance between the spins. 
Interest in the RVB state revived after the discovery of high temperature superconductivity in 1987[^. 
The common structural ingredient of the high— Tc cuprate systems is the copper-oxide {CUO2) plane which 
ideally behaves as a S" = ^ HAFM defined on a square lattice. It is largely agreed that the ground state 
(T = 0) has AFM long range order (LRO). The low-lying excitations are the conventional S = 1 magnons. 
In the spinon picture, a magnon is a pair of confined spinons. The spinons cannot move apart from each 
other unlike in Id. The cuprates exhibit a rich phase diagram as a function of the dopant concentration. On 
doping, positively charged holes are introduced in the Cu02 plane. The holes are mobile in a background 
of antiferromagnetically interacting spins. The motion of holes acts against antiferromagnetism and the 
AFM LRO is rapidly destroyed as the concentration of holes increases. The resulting spin disordered state 
has been speculated to be a RVB state. In close analogy with the S — ^ HAFM chain, the low-lying spin 
excitations in the RVB state are pairs of spinons. The spinons are created by breaking a VB. The spinons 
are not confined as in the case of an ordered ground state but separate via a rearrangement of the VBs. The 
spinons have spin ^ and charge 0. The charge excitations in a RVB state are known as holons with charge 
+e and spin 0. Holons are created on doping the RVB state, i.e., replacing electrons by holes. Spinons and 
holons are best decribed as topological excitations in a QSL. The key feature of the doped RVB state is that 
of spin-charge separation, i.e., the spin and charge excitations are decoupled entities. Spin-charge separation 
can be rigorously demonstrated in the case of interacting electron systems in Id known by the general name 
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of Luttinger Liquids (LLs). The Hubbard model in Id is the most well-known example of a LL. The model 
is integrable and the BA results for the excitation spectrum confirm that the spinous and the holons are the 
elementary excitations ||36| , |37| . 

Coming back to the RVB state, there has been an intensive search for spin models in 2d with RVB states 
as exact ground states. Recent calculations show that there is AFM LRO in the ground state of the S — ^ 
HAFM on the triangular lattice, contrary to Anderson's original conjecture ||3^. Frustrated spin models with 
n.n. as well as non-n.n. exchange interactions have been constructed for which the RVB states are the exact 
ground states in certain parameter regimes |Q. These are short-ranged RVB states with the VBs forming 
between n.n. spin pairs. The spinon excitation spectrum in this case is gapped. A model which captures 
the low energy dynamics in the RVB scenario is the Quantum Dimer Model (QDM)||40). The Hamiltonian 
of the model defined on a square lattice is given by 



H„oM = x{-t(|i ^(::|. H.c) 
-Kl::>(::l Hi iXi il)) 



(23) 



where the solid lines represent dimers (VBs) and the sum runs over all the plaquettes of the lattice. 
The first term of the Hamiltonian is the kinetic part representing the fiipping of a pair of parallel dimers 
on the two bonds of a plaquette to the other possible orientation, i.e., from horizontal to vertical and vice 
versa. The second term counts the number of fiippable pairs of dimers in any dimer configuration and is 
analogous to the potential term of the Hamiltonian. The ground state of the QDM on the square lattice 
is not, however, a QSL except at the special point t — V. Moessner and Sondhi ||4l) have studied the 
QDM on the triangular lattice and shown that, in contrast to the square lattice case , the ground state is 
a RVB state with deconfined, gapped spinons in a finite range of parameters. Recently, some microscopic 
models of 2d magnets have been proposed |Q the low-lying excitations of which are of three types: spinons, 
holons and "vortex-like" excitations with no spin and charge, dubbed as visons. Some of these models are 
related to the QDM. Two integrable models |43| have been constructed which share common topological 
features with the microscopic models in 2d and have applications in fault-tolerant quantum computation. 
The models , however, cannot resolve the issue of spinons in 2d as quantum numbers like the total are not 
conserved in these models. The search for microscopic models in 2d, with spinons as elementary excitations, 
acquires particular significance in the light of recent experimental evidence of the spinon continuum in the 
2d frustrated quantum antiferromagnet Cs2CuCli |Q. The ground state of this compound is expected to 
be a QSL with spinons and not magnons as elementary excitations. Exactly solvable models in 2d are needed 
for a clear understanding of the origin of the experimentally observed spinon continuum. 

Real materials are often anisotropic in character. The anisotropy may be present in the exchange interac- 
tion Hamiltonian itself or there may be additional terms in the Hamiltonian corresponding to different types 
of anisotropy. A well-known anisotropic interaction, present in many AFM materials, is the Dzyaloshinskii- 
Moriya (DM) interaction with the general form 

Hdm = X ^j) (24) 

Moriya ji^ provided the microscopic basis of the DM interaction by extending Anderson's superexchange 
theory to include the spin-orbit interaction. The DM coupling acts to cant the spins because the coupling 
energy is minimised when the two spins are perpendicular to each other. Some examples of materials with DM 
interaction include the quasi-2d compound Cs2CuCl i p^ , the Cu02 planes of the undoped cuprate system 
La2CuOi the quasi-ld compound Cu-Benzoate ||47[| etc. The DM canting of spins is responsible for the 
small ferromagnetic moment of the Cu02 planes even though the dominant in-plane exchange interaction is 
AFM in nature. Alcaraz and Wreszinski have shown that the XXZ quantum Heisenberg chain (both FM 
and AFM) with DM interaction is equivalent to the XXZ Hamiltonian with modified boundary conditions 
and anisotropy parameter 4^ . The DM interaction is assumed to be of the form 
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A ^ 

Hdm{^) = -- ^«af+^ _ <^ J (25) 

1=1 

,i.e., the vector l5 in Eq.(23) is in the z-direction. The new anisotropy parameter is ^^^^^ where 5 is the 
anisotropy parameter of the original XX Z Hamiltonian. With changed boundary conditions, the model is 
still BA solvable. In fact, in the thermodynamic limit (iV oo), the boundary conditions do not affect 
the critical behaviour. Thus, the Hamiltonian, which includes both the XXZ Hamiltonian and the DM 
interaction, has the same critical properties and the phase diagram as the XXZ Hamiltonian with the 
anisotropy parameter --^=^. 

We next turn our attention to spin— S" {S > ^) quantum spin chains. The spin— S* Heisenberg exchange 
interaction Hamiltonian in Id is not integrable. A family of Heisenberg-like models has been constructed for 
S — 1,^, 2, |, ... etc. for which the spin— S* quantum Hamiltonian is given by 

=^Q(^,.^,+i) (26) 

i 

where Q{x) is a polynomial of degree 25* |^^. With this generalization, the spin— S* quantum spin chains 
are integrable. The integrable models, however, do not distinguish between half-odd integer and integer 
spins. In both the cases, the integrable models have gapless excitation spectrum. For half-odd integer 
AFM Heisenber g sp in chains (with only the bilinear exchange interaction term), the Lieb-Schultz-Mattis 
(LSM) theorem ||50|| states that the excitation spectrum is gapless. The theorem cannot be proved for AFM 
integer spin chains. Haldane in 1983 pointed out the difference between the half-odd integer and integer AFM 
Heisenberg spin chains and made the conjecture that integer spin chains have a gap in the excitation spectrum 
Integer spin quantum antiferromagnets in Id have been widely studied analytically, numerically and 
experimentally and Haldane's conjecture has turned out to be true. There are several examples of quasi- 
Id S* = 1 AFM materials which exhibit the Haldane gap. Some of the most widely studied materials are 
CsNiCh, Ni{C2HsN2)2N02{ClOi) (NENP), Y2BaNiO^ etc. Recently, experimental evidence of a S" = 2 
antiferromagnet which exhibits the Haldane gap has been obtained. In this compound the manganese ions 
form effective S = 2 spins and are coupled in a quasi-ld chain ||5^. Integrable models of integer spin chains 
do not reproduce the Haldane gap but are of considerable interest since they provide exact information about 
the phase diagram of generalised integer spin models. Consider the generalised Hamiltonian for an AFM 
S = \ chain: 



H = Y, [cose{^ ,r§ + 5zn0( 5^,. S^,+i)2] (27) 



with 9 varying between and 27r. The biquadratic term has been found to be relevant in some real integer- 
spin materials. There are two gapped phases: the Haldane phase for — ^ < 6* < |-and a dimerised phase for 
— ^ < 6* < — ^ [Q. At 9 = the model is integrable and the gap vanishes to zero. This point separates 
the two gapped phases, Haldane and dimerised, which have different symmetry properties. Thus a quantum 
phase transition occurs aX d — from the Haldane to the dimerised phase. The integrable model provides 
exact location of the transition point. The point 9 — ^ corresponds to the Hamiltonian which is a sum over 
permutation operators and is again exactly solvable. The Haldane phase includes the isotropic Heisenberg 
chain {9 = 0) and the AfHeck-Kennedy-Lieb-Tasaki (AKLT) Hamiltonian {tan9vBS = |) The latter 

model is not integrable but the ground state is known exactly. The ground state is described as a valence 
bond solid (VBS) state in which a VB (singlet) covers every link of the chain. Since the gap does not become 
zero for < 9 < 9vbs, there is no phase transition in going from one limiting Hamiltonian to the other. 
Thus the isotropic Heisenberg and AKLT chains are in the same phase. 

The doped cuprate systems exhibit a variety of novel phenomena in their insulating, metaUic and su- 
perconducting phases. A full understanding of these phenomena is as yet lacking. There is currently a 
strong research interest in doped spin systems. The idea is to look for simpler spin systems in which the 
consequences of doping can be studied in a less ambiguous manner. The spin— 1 HG nickelate compound 
Y2BaNiOz can be doped with holes on replacing the off-chain ions by Ca^+ ions. Inelastic neutron 
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scatterin g ( INS) measurements on the doped compound provide evidence for the appearance of new states in 
the HG 1^^. The structure factor S{q), obtained by integrating the dynamical correlation function S{q,uj) 
over oj, acquires an incommensurate, double-peaked form in the doped state |^^. Frahm et al have con- 
structed an integrable model describing a doped spin— 1 chain. In the undoped Hmit, the spectrum is gapless 
and so the HG of the integer spin system is not reproduced. It is, however, possible to reintroduce a gap 
in the continuum limit where a field-theoretical description of the model is possible. The model has limited 
relevance in explaining the physical features of the doped nickelate compound. Another interesting study 
relates to the appearance of magnetization plateaus in the doped S = 1 integrable model . The location 
of the plateaus depends on the concentration of holes. Experimental evidence of this novel phenomenon has 
not been obtained so far. 

An electron in a sohd, locahsed around an atomic site, has three degrees of freedom charge, spin and 
orbital. The orbital degree of freedom is relevant to several transition metal oxides which include the 
cuprate and manganite systems. The latter compounds on doping exhibit the phenomenon of colossal 
magnetoresistance in which there is a huge change in electrical resistivity on the application of a magnetic 
field. The manganites like the cuprates have a rich phase diagram as a function of the dopant concentration 
| p9| . We now give a specific example of the orbital degree of freedom. The Mn^+ ion in the manganite 
compound LaMnO^ has four electrons in the outermost 3d energy level. The electrostatic field of the 
neighbouring oxygen ions splits the 3d energy level into two sublevels, t2g and Cg. Three of the four electrons 
occupy the three t2g orbitals dxy,dyzdzx and the fourth electron goes to the e^-sublevel containing the 
two orbitals d^-^-y-i and d^^^-r'^. The fourth electron thus has an orbital degree of freedom as it has two 
possible choices for occupying an orbital. The four electrons have the same spin orientation to minimise the 
electrostatic repulsion energy according to the Hund's rule. The total spin is thus S — 2. The orbital degree 
of freedom is described by the pseudospin 7^ such that = \{—\) when the dx^_y2 {d^z^_^2) orbital is 
occupied. The three components of the pseudospin satisfy commutation relations similar to those of the spin 
components. The Cg doublet is further split into two hyperfine energy levels due to the well-known Jahn- Teller 
(JT) effect. In concentrated systems, the JT effect can lead to orbital ordering below an ordering temperature. 
In the antiferromagnetically ordered Neel state, the spins are alternately up and down. Similarly, in the case 
of antiferroorbital ordering, the occupied orbitals alternate in type at successive sites of the lattice. The 
orbital degree of freedom is frozen as a result. Apart from the JT mechanism of orbital ordering, there is 
an exchange mechanism which may lead to orbital order. The exchange mechanism is a generalisation of 
the usual superexchange to the case of orbital degeneracy. Starting from the degenerate Hubbard model, 
in which there are two degenerate orbitals at each site, one can derive the following generalised exchange 
Hamiltonian | |60| : 

H = Y, {Jl~^^^l + J2^^^l + J3(5^.. 5^,)(T^..T^,)} (28) 

Consider the case Ji — J2 — J- For J3 = 0, two independent Heisenberg-like Hamiltonians are obtained 
which are BA solvable. At the Kolezhuk-Mikeska point, ^ = |, the ground state is exactly known |Q. The 
point ^ = 4 is integrable and there are three gapless excitation modes. The compounds Na2Ti2Sb20 and 
NaV205 are examples of materials in Id with coupled spin and orbital degrees of freedom ||6^. These systems 
have been described by anisotropic versions of the Hamiltonian in Eq.(28) but without adequate agreement 
with experiments. The elementary excitations in the orbital sector are the orbital waves or "orbitons". An 
excitation of this type is created in the orbitally ordered state by changing the occupied orbital at a site 
and letting the defect propagate in the solid. The excitations are analogous to the spin waves or magnons 
in a magnetically ordered solid. Experimental evidence of orbital waves has recently been obtained in the 
manganite compound LaMnO^ through Raman scattering measurements | |6^ , |64[ . As discussed before, 
integrable spin models provide important links between theory and experiments. A similar scenario in the 
case of systems with coupled spin and orbital degrees of freedom is yet to develop. 
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3 Ladder models 
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Figure 2. ^4 two-chain ladder. The rung and intra-chain n.n. exchange interactions 
are of strength Jr and J respectively. 



The simplest ladder model consists of two chains coupled by rungs (Figure 2). In general, the ladder 
may consist of n chains coupled by rungs. In the spin ladder model, each site of the ladder is occupied by 
a spin (in general of magnitude ^) and the spins interact via the Heisenberg AFM exchange interaction. 
In the doped spin ladder model, some of the sites are empty, i.e., occupied by holes. The holes can move 
in the background of interacting spins. There are two major reasons for the considerable research interest 
in ladders. Powerful techniques like the BA and bosonization are available for the study of Id many body 
systems whereas practically very few rigorous results are known for 2d systems. Ladders provide a bridge 
between Id and 2d physics and are ideally suited to study how the electronic and magnetic properties change 
as one goes from a single chain to the square lattice. The unconventional properties of the Cu02 planes of 
the cuprate systems are the main reason for the significant interest in 2d many body systems. Many of these 
properties are ascribed to strong correlation effects. Ladders are simpler systems in which some of the issues 
associated with strong correlation can be addressed in a more rigorous manner. The second motivation for 
the study of ladder systems is that several such systems have been discovered in the recent past. In the 
following, we describe in brief some of the major physical properties of ladders. There are two exhaustive 
reviews on ladders which provide more detailed information ||65| , |66[ . 

Consider a two-chain spin ladder described by the AFM Heisenberg exchange interaction Hamiltonian 

H = Y,J^i^^^l (29) 

The n.n. intra-chain and the rung exchange interactions are of strength J and Jr respectively. When Jr — 0, 
one obtains two decoupled AFM spin chains for which the excitation spectrum is known to be gapless. For 
all ^ > 0, a gap (the so-called spin gap (SG)) opens up in the spin excitation spectrum. The result is 
easy to understand in the simple limit in which the exchange coupling Jr along the rungs is much stronger 
than the coupHng J along the chains. The intra-chain coupHng may thus be treated as perturbation. When 
J = 0, the exact ground state consists of singlets along the rungs. The ground state energy is —^^f^, where 
N is the number of rungs in the ladder. The ground state has total spin 5 = 0. In first order perturbation 
theory, the correction to the ground state energy is zero. A S" = 1 excitation may be created by promoting 
one of the rung singlets to a 5' = 1 triplet. The weak coupling along the chains gives rise to a propagating 
S = 1 magnon. In first order perturbation theory, the dispersion relation is 

t^(fc) = Jn. + Jcosk (30) 
where k is the momentum wave vector. The SG defined as the minimum excitation energy is given by 



^SG 



uj{tt) ~ {Jr - J) (31) 



The two-spin correlations decay exponentially along the chains showing that the ground state is a quantum 
spin liquid (QSL). The magnons can further form bound states. Experimental evidence of two-magnon 
bound states has been obtained in the S = ^two-chain ladder compound Cai4-xLaxCu2404i {x = 5 and 4) 
|R7|. The family of compounds Srn-iCun+i02n consists of planes of weakly-coupled ladders of ^^chains 
||68 |. For n = 3 and 5, respectively, one gets the two-chain and three-chain ladder compounds SrCuiO^ and 
Sr2Cu^0z respectively. For the first compound, experimental evidence of the SG has been obtained. The 
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latter compound has properties similar to those of the Id Heisenberg AFM chain |^^. A recent example of a 
spin ladder belonging to the organic family of materials is the compound {C5Hi2N)2CuBr4, a ladder system 
with strong rung coupling ~ 3.5) |^^. The phase diagram of the AFM spin ladder in the presence of 
an external magnetic field is particularly interesting. In the absence of the magnetic field and at T = 0, 
the ground state is a QSL with a gap in the excitation spectrum. At a field i?ci , there is a transition to 
a gapless Luttinger Liquid (LL) phase (g/isiJci = ^SG, the spin gap, fiB is the Bohr magneton and g the 
Lande splitting factor). There is another transition at an upper critical field to a fully polarised FM 
state. Both Hc^ and i?c2 are quantum critical points. The quantum phase transition from one ground state 
to another is brought about by changing the magnetic field. At small temperatures, the behaviour of the 
system is determined by the crossover between two types of critical behaviour: quantum critical behaviour 
at T = and classical critical behaviour at T 0. Quantum effects are persistent in the crossover region at 
small finite temperatures and such effects can be probed experimentally. In the case of the ladder system 
{C5Hi2N)2CuBr4, the magnetization data obtained experimentally exhibit universal scaling behaviour in 
the vicinity of the critical fields i/ci and ■ In the gapless regime i?ci < H < , the ladder model can be 
mapped onto an XX Z chain the thermodynamic properties of which can be calculated exactly by the BA. 
The theoretically computed magnetization M versus magnetic field h curve is in excellent agreement with 
the experimental data. Organic spin ladders provide ideal testing grounds for the theories of quantum phase 
transitions. For inorganic spin ladder systems, the value of ifci is too high to be experimentally accessible. 

Bose and Gayen ||7^ have studied a frustrated two-chain spin model with diagonal couplings. The 
intrachain and diagonal spin-spin interactions are of equal strength J. It is easy to show that for > 2 J, 
the exact ground state consists of singlets (dimers) along the rungs with the energy Eg = — MaiY. where TV is 
the number of rungs. Xian j?^ later pointed out that as long as ^ > {^)^ — 1-401, the rung dimer state 
is the exact ground state. At ^ = (4f) ' there is a first order transition from the rung dimer state to the 
Haldane phase of the S — I chain. Kolezhuk and Mikeska |Q have constructed a class of generalised S = ^ 
two-chain ladder models for which the ground state can be determined exactly. The Hamiltonian H is a 
sum over plaquette Hamiltonians and each such Hamiltonian contains various two-spin as well as four-spin 
interaction terms. They have further introduced a toy model which has a rich phase diagram in which the 
phase boundaries can be determined exactly. 

The standard spin ladder models with bilinear exchange are not integrable. For integrability, multispin 
interaction terms have to be included in the Hamiltonian. Some integrable ladder models have already been 
constructed [Q. We discuss one particular model proposed by Wang ||7^. The Hamiltonian is given by 



Ji ^ J2 ^ 

JJ ^ JT ^ 

+ T ^ (^^■•^^■+1) (^.■■^.+1) + f ^ (32) 

where and j are the PauH matrices associated with the site j of the upper and lower chains respectively. 
A'' is the total number of rungs in the system. The ordinary spin ladder Hamiltonian is obtained from Eq. 
(32) when the four spin terms are absent, i.e., Ui = U2 = 0. For general parameters Ji, J2, Ui and U2, the 
model is non-integrable. The integrable cases correspond to C/i = Ji, t/2 = or Ui = Ji,U2 — Without 
loss of generality one can put Ji = Ui = 1, J2 = J and U2 = U. For U the Hamiltonian (32) reduces to 

1 ^ 7 ^ 11 

^ = i + ^3-^3 + l)(^ + ^3--^ 3 + 1) + 2 E(^^-^^- - 1) + - (33) 

Three quantum phases are possible. For J > J'^ = 2, the system exists in the rung dimerised phase. The 
ground state is a product of singlet rungs. The SG is given by Asa = 2(J — 2). For > J > JI, a 
gapless phase is obtained with three branches of gapless excitations. is the quantum critical point at 
which a QPT from the dimerised phase to the gapless phase occurs. In the vicinity of the quantum critical 
point, the susceptibility and the specific heat can be calculated using the thermodynamic BA. From the 
low-temperature expansion of the thermodynamic BA equation, one obtains 
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(34) 



which are typical of quantum critical behaviour. In the presence of an external magnetic field h, the magnetic 
field can be tuned to drive a QPT at the quantum critical point he — 2(J — 2) from the gapless phase to a 
gapped phase. The third quantum phase {h = 0) is obtained for J < Jl — + This is a gapless 

phase with two branches of gapless excitations. For U = — i, a similar phase diagram is obtained. Note 
that the ladder model may equivalently be considered as a spin-orbital model with V and representing 
the spin and the pseudospin. 

Doped ladder models are toy models of strongly correlated systems |Q. In these systems, the double 
occupancy of a site by two electrons, one with spin up and the other with spin down, is prohibited due to 
strong coulomb correlations. In a doped spin system, there is a competition between two processes: hole 
delocalization and exchange energy minimization. A hole moving in an antiferromagnetically ordered spin 
background, say, the Neel state, gives rise to parallel spin pairs which raise the exchange interaction energy 
of the system. The questions of interest are: whether a coherent motion of the holes is possible, whether 
two holes can form a bound state, the development of superconducting (SC) correlations, the possibility 
of phase separation of holes etc. Some of these issues are of significant relevance in the context of doped 
cuprate systems in which charge transport occurs through the motion of holes |Q. In the SC phase, the 
holes form bound pairs with possibly d-wave symmetry. Several proposals have been made so far on the 
origin of hole binding but there is as yet no general consensus on the actual binding mechanism. The doped 
cuprate systems exist in a 'pseudogap' phase before the SC phase is entered. In fact, some cuprate systems 
also exhibit SG. As already mentioned, the doped two-chain ladder systems are characterised by a SG. The 
issue of how the gap evolves on doping is of significant interest. The possibility of binding of hole pairs in a 
two-chain ladder system was first pointed out by Dagotto et al In this case, the binding mechanism is 
not controversial and can be understood in a simple physical picture. Again, consider the case Jr ^ J, i.e., 
a ladder with dominant exchange interactions along the rungs. In the ground state, the rungs are mostly 
in singlet spin configurations. On the introduction of a single hole, a singlet spin pair is broken and the 
corresponding exchange interaction energy is lost. When two holes are present, they prefer to be on the 
same rung to minimise the loss in the exchange interaction energy. The holes thus form a bound pair. In 
the more general case, detailed energy considerations show that the two holes tend to be close to each other 
and effectively form a bound pair. For more than two holes, several calculations suggest that considerable 
SC pairing correlations develop in the system on doping. True superconductivity can be obtained only in 
the bulk limit. Theoretical predictions motivated the search for ladder compounds which can be doped with 
holes. Much excitement was created in 1996 when the ladder compound S'ri4_a;Caa;Cu2404i was found to 
become SC under pressure at x — 13.6 j?^. The transition temperature Tc is ^ 12K at a pressure of 3GPa. 
As in the case of cuprate systems, bound pairs of holes are responsible for charge transport in the SC phase. 
Experimental results on doped ladder compounds point out strong analogies between the doped ladder and 
cuprate systems []65| . 

The strongly correlated doped ladder system is described by the t-J Hamiltonian 

Ht-j = -J2 t^d^^A + H.C.) + J,, (^..^, - -n^Uj) (35) 

The and Ci^ are the electron creation and annihilation operators which act in the reduced Hilbert space 
(no double occupancy of sites), 

5+ = C+(l-n,_,) 

a, = Q„(l-n,_,) (36) 

where cr is the spin index and n^, nj are the occupation numbers of the ith and jth sites respectively. The 
first term in Eq.(35) describes the motion of holes with hopping integrals and t for motion along the 
rung and chain respectively. In the standard t — J ladder model, i and j are n.n. sites. The second term 
contains the usual AFM Heisenberg exchange interaction Hamiltonian. The t — J model thus describes the 
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motion of holes in a background of antiferromagnetically interacting spins. A large number of studies have 
been carried out on t — J ladder models. These are reviewed in Refs. ||6^ , |6^ . We describe briefly some of 
the major results. The SG of the undoped ladder changes discontinuously on doping. Remember that the 
SG is the difference in energies of the lowest triplet excitation and the ground state. In the doped state, 
there are two distinct triplet excitations. One triplet excitation is that of the undoped ladder obtained by 
exciting a rung singlet to a rung triplet. A new type of triplet excitation is possible when at least two holes 
are present. On the introduction of two holes in two rung singlets, a pair of free spin— ^'s is obtained which 
combines to give rise to a singlet {S = 0) or a triplet {S = 1) state. The triplet conflguration of the two free 
spins corresponds to the second type of triplet excitation. The SG of this new excitation is unrelated to the 
SG of the magnon excitation. The true SG is the one which has the lowest value in a particular parameter 
regime. 

The low energy modes of a ladder system are characterised by their spin. Singlet and triplet excitations 
correspond to charge and spin modes respectively. In each sector, the hole may further be in a bonding 
or antibonding state with opposite parities. We consider only the even parity sector to which the lowest 
energy excitations belong. In both the S = and S — 1 sectors, an excitation continuum with well-defined 
boundaries is present. The S = and the S = 1 continua are degnerate in energy. A bound state branch 
with S = splits off below the continuum the lowest energy of which corresponds to the c.o.m. momentum 



wave vector K = |79, Thus the two-hole ground state is in the singlet sector and corresponds to a 
bound state of two holes with K = 0. The bound state has d— wave type symmetry. Within the bound state 
branch, excitations with energy infinitesimally close to the ground state are possible. These excitations are 
the charge excitations since the total spin is still zero and the charge excitation spectrum is gapless. The 
lowest spin excitations in a wide parameter regime are between the S — ground state and the lowest energy 
state in the S — 1 continuum js^. The continuum does not exist in the undoped ladder and so the SG 
evolves discontinuously on doping in this parameter regime. A suggestion has, however, been made that the 
lowest triplet excitation is a bound state of a magnon with a pair of holes |^. In summary, the two-chain 
ladder model has the feature that the charge excitation is gapless but the spin excitation has a gap. This is 
the Luther-Emery phase and is different from the LL phase in which both the spin and charge excitations 
are gapless. 

Bose and Gayen have derived several exact, analytical results for the ground state energy and the low- 
lying excitation spectrum of the frustrated t— J ladder doped with one and two holes. The undoped ladder 
model has already been described. In the doped case, the hopping integral has the value tR for hole motion 
along the rungs and the intra-chain and diagonal hopping integrals are of equal strength t. The latter 
assumption is crucial for the exact solvability of the eigenvalue problem in the one and two hole sectors. 
Though the model differs from the standard t — J ladder model (the diagonal couplings are missing in the 
latter), the spin and charge excitation spectra exhibit similar features. In particular, the dispersion relation 
of the two-hole bound state branch is obtained exactly and the exact ground state is shown to be a bound 
state of two holes with K — and d— wave type symmetry. The ladder exists in the Luther-Emery phase. 
There is no spin charge separation, as in the case of a LL. In the exact hole eigenstates, the hole is always 
accompanied by a free spin— i. The hole-hole correlation function can also be calculated exactly. When 
Jr ^ J, the holes of a bound pair are predominantly on the same rung. For lower values of Jr, the holes 
prefer to be on n.n. rungs so that energy gain through the delocahzation of a hole along the rung is possible. 

The t-J ladder model constructed by Bose and Gayen is not integrable. Frahm and Kundu jsj] have 
constructed at — J ladder model which is integrable. The Hamiltonian is given by 

H = J2 Ht-\j + H^nt + Hrung - fiU (37) 
a 

The two chains of the ladder are labelled by a = 1,2 and /i is the chemical potential couphng to the number of 
electrons in the system. Hj:'^j is the t— J Hamiltonian (Eq.(35)) for a chain plus the terms n^"' -f nj'^^-^ where 
n^-^'is the total number of electrons on site j. 



Hi, 



(38) 
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Hrung includes the t — J Hamiltonian (Eq. (35)) corresponding to a rung and a Coulomb interaction term 
VJ2j '^i^^'^j^^- The possible basis states of a rung are the following. When no hole is present, a rung can be 
in a singlet or a triplet spin configuration. When a single hole is present, the rung is in a bonding (|cr+)) or 
antibonding {\<7-)) state with \a±) = (WO) ± |Ocr)) and cr =1 or |. The rung can further be occupied by 
two holes. Frahm and Kundu have studied the phase diagram of the ladder model at low temperatures and 
in the strong coupling regime Jr :s> l,V + \tji\ near half-filling. In this regime, the triplet states are 
unfavourable. By excluding the triplet states and choosing J = 2t ~ 2, the Hamiltonian H (Eq.(37)) can be 
rewritten as 

5 

H = -^I{yj+i-^AiNi-\' const. (39) 
j 1=1 

where Ni,l = 1, 2(3, 4) is the number of bonding (antibonding) single hole rung states with spin f , i and TVs 
is the number of empty rungs. If L is the total number of rungs in the ladder, the remaining Nq — L — J2i 
rungs are in singlet spin configurations. The permutation operator Ujk interchanges the states on rungs j 
and k. If both the rungs are singly occupied by a hole, an additional minus sign is obtained on interchanging 
the rung states. The potentials j4;'s are: 



Ai=A2 = n+=tB.-H + V (40) 

A3=Ai = n^=-tR^fi + V (41) 

A5 = V ^ -2n+V (42) 

The nature of the ground state and the low-lying excitation spectrum depends on the relative strengths of 
the potentials Ai's. The Hamiltonian (39) is BA solvable. The phase diagram V vs. the hole concentration 
Uh has been computed for /i+ = ^_,i.e., tu — 0. For large repulsive V, the ground state can be described 
as a Fermi sea of single hole states \a±) propagating in a background of rung dimer states \s). The double- 
hole rung states \d) are energetically favourable for sufficiently strong attractive rung interactions. In the 
intermediate region, both types of hole rung states are present. In the frustrated t — J ladder model studied 
by Bose and Gayen |Q , the exact two-hole ground state is a linear combination of single-hole and double-hole 
rung states propagating in a background of rung dimer states. The single-hole rung states are the bonding 
states. 

In a remarkable paper, Lin et al have considered the problem of electrons hopping on a two-chain 
ladder. The interaction between the electrons is sufficiently weak and finite-ranged. At half filling, a 
perturbative renormalization group (RG) calculation shows that the model scales onto the Gross-Neveu 
(GN) model which is integrable and has SO{8) symmetry. At half filling, the two-chain ladder is in the Mott 
insulating phase with d-wave pairing correlations. The insulating phase is further a QSL. The integrability 
has been utilised to determine the exact energies and quantum numbers of all the low energy excitations 
which constitute the degenerate SO{8) multiplets. The lowest-lying excitations can be divided into three 
octets all with a non-zero gap (mass gap) m. Each excitation has a dispersion ei{q) = \Jm?- + where q 
is the momentum variable measured w.r.t. the minimum energy value. One octet consists of two-particle 
excitations: two charge ±2e Cooper pairs around zero momentum, a triplet of S* = 1 magnons around 
momentum (tt, tt) and three neutral S" = particle-hole pair excitations. 50(8) transformations rotate the 
components of the vector multiplet into one another unifying the excitations in the process. The 5*0(5) 
subgroup which rotates only the first five components of the vector is the symmetry proposed by Zhang |Q 
to unify antiferromagnetism and superconductivity in the cuprates. The vector octet is related by a triality 
symmetry to two other octets with mass gap m. The 16 particles of these two octets have the features of 
quasi-electrons and quasi-holes. Above the 24 states with mass gap m, there are other higher- lying "bound" 
states with mass gap VSm. Finally, the continuum of scattering states occurs above the energy 2ra. Lin et 
al has further studied the effects of doping a small concentration of holes into the Mott insulating phase. 
In this Hmit, the effect of doping can be incorporated in the GN model by adding a term —^iQ to the 
Hamiltonian, /i being the chemical potential and Q the total charge. Integrability of the GN model is not 
lost as Q is a global 50(8) generator. Doping is possible only for 2/i > m when Cooper pairs enter the 
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system. The doped ladder exists in the Luther-Emery phase, whereas in the half-filled insulating limit both 
the spin and charge excitations are gapped. In the doped phase, the Cooper pairs can transport charge and 
quasi-long-range d— wave SC pairing correlations develop in the system. The other features of the standard 
t— J ladder model, e.g., the discontinuous evolution of the SG on doping is reproduced. The lowest triplet 
excitation is a bound state of a = 1 magnon with a Cooper pair. As mentioned before, a similar result has 
been obtained numerically in the case of the standard t — J ladder |^2|. The triplet excitation belongs to 
the family of 28 excitations with mass gap V3™. If x denotes the dopant concentration, then the SG jumps 
from As(a; 0) = to to As(a; = 0+) = (^3-1) m upon doping. The integrability of the weakly-interacting 
two-chain ladder model has yielded a plethora of exact results which illustrate the rich physics associated 
with undoped and doped ladders. 

4 Concluding Remarks 

Integrable models have a dual utility. They serve as testing grounds for approximate methods and techniques. 
Also, they are often models of real systems and provide rigorous information about the physical properties 
of such systems. Integrable models are sometimes more general than what are required to describe real 
systems. In such cases, an integrable model corresponds to an exactly solvable point in the general phase 
diagram. The point may be a quantum critical point at which transition from one quantum phase to another 
occurs or the integrable model may be in the same phase as a more realistic model. In the latter case, the 
physical properties of the two models are similar. In this review, we have discussed the physical basis of some 
integrable spin models with special focus on the relevance of the models to real systems. The Heisenberg spin 
chain is probably the best example of the essential role played by exact solvability in correctly interpreting 
the experimental data. The concept of spinous owes its origin to the exact analysis of the BA equations. The 
theoretical prediction motivated the search for real spin systems in which experimental confirmation could 
be made. In this review, examples are also given of systems for which the links between integrable models 
and experimental results are not well established. A major portion of the review is devoted to physical 
systems which exhibit rich phenomena, like the systems with both spin and orbital degrees of freedom 
and undoped and doped spin ladder systems, where the need for integrable systems is particularly strong. 
These systems exhibit a variety of novel phenomena a proper understanding of which should be based on 
rigorous theory. Two-dimensional spin systems with QSL ground states have been specially mentioned to 
explain the recent interest in constructing integrable models of such systems. The review is meant to be an 
elementary introduction to the genesis and usefulness of integrable models vis-a-vis physical spin systems. 
Future challenges are also highlighted to motivate further research on integrable models. 

There are some AFM spin models which are not integrable but for which the ground states and in some 
cases the low-lying excited states are known exactly. The most prominent amongst these are the Majumdar- 
Ghosh (MG) chain and the AKLT jsij model respectively. The MG Hamiltonian is defined in Id for 
spins of magnitude \. The Hamiltonian includes both n.n. as well as n.n.n. interactions. The strength of the 
latter is half that of the former. The exact ground state is doubly degenerate and the states consist of singlets 
along alternate links of the lattice. The excitation spectrum is not exactly known and has been calculated 
on the basis of a variational wave function ||8^ . Generalizations of the MG model to 2d with exactly-known 
ground states are possible |^, |9[ |9^, 0. The Shastry-Sutherland model |^ is of much current interest 
due to the recent discovery of the compound SrCu2{BO^)2 which is well-described by the model [Q. Some 
of these models including the AKLT model have been reviewed in the references fo^ , |9^ from which 

more information about the models can be obtained. These models incorporate physical features of real 
systems and provide valuable insight on the magnetic properties of low-dimensional quantum spin systems. 
The models supplement integrable models in obtaining exact information and provide motivation for the 
construction of integrable generalisations. 
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